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Second order brane cosmology with radion stabilization 
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We study cosmology in the five-dimensional Randall-Sundrum brane-world with a stabilizing 
effective potential for the radion and matter localized on the branes. The analysis is performed 
by employing a perturbative expansion in the ratio p/V between the matter energy density on the 
branes and the brane tensions around the static Randall-Sundrum solution (which has p = 0 and 
brane tensions ±R). This approach ensures that the matter evolves adiabatically and allows us 
to find approximate solutions to second order in p/V. Some particular cases are then analyzed in 
details. 


PACS numbers: 12.60.-i, 04.50.+h, 98.80.Cq, 98.80.-k 

I. INTRODUCTION 

Higher and higher precision data which are about to 
be collected in new experiments of particle physics and 
astrophysics in the next few years convey considerable 
attention to theories with extra dimensions. The main 
role of such theories, originally introduced in the 20’s 
by Kaluza and Klein P, is to provide a connection 
between particle physics and gravity at some level. At 
a deeper level, string theory unifies all the interactions 
by means of some n-dimensional manifold (with n > 4) 
where the fundamental objects are supposedly living; at 
a more phenomenological level, models which assume the 
existence of extra dimensions, no matter their origin, 
are considered in order to solve some puzzles of parti¬ 
cle physics, cosmology and astrophysics, giving rise to 
many possible observable consequences. 

Originally proposed in order to solve the problem of 
the large hierarchy between Gravity and Standard Model 
scales, the Randall-Sundrum model of Ref. 0 (RS I) has 
acquired considerable relevance due to its stringy inspira¬ 
tion. It represents the prototype of the so-called brane- 
world and differs from previous models in that it con¬ 
strains standard matter on a four-dimensional manifold 
(the brane) just letting gravity (and exotic matter) prop¬ 
agate everywhere. The RS I solution to the hierarchy 
problem needs one additional compactified (orbifolded) 
spatial dimension with two branes located at its fixed 
points, plus a negative cosmological constant filling the 
space between such branes (the bulk). The bulk cosmo¬ 
logical constant A warps the extra dimension and gen¬ 
erates the effective four-dimensional physical constants 
we measure. It was soon realized that the modifications 
to four-dimensional gravity induced by the fifth dimen¬ 
sion may be reduced to such a short distance effect to 
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be unobservable even in the presence of just one brane 
and infinite compactification radius (the RS II model of 

Ref. 0 ). 

The cosmological features of the RS models are nowa¬ 
days being investigated even more than its particle 
physics consequences, due to the refined results lately 
obtained and to the major problems recent astrophysi- 
cal data have revealed: the possible late time accelera¬ 
tion from supernovae, GMBR spectrum, dark matter and 
dark energy quests suggest either a full revision of the 
modern theoretical physics approach or the possibility of 
the existence of further, up to now ignored, ingredients 
such as the extra dimensions. 

In particular the single brane RS II cosmological dy¬ 
namics Hi is known to generate {p/V)^ corrections to 
standard Friedmann and acceleration equations, where p 
is the energy density of the fluid filling the brane and V 
is the constant brane tension. These corrections are neg¬ 
ligible when p <^V, the regime in which the RS II model 
is reliable and leads to standard cosmic evolution. The 
two brane RS I setup is much more involved: a stabi¬ 
lization mechanism for the distance between the branes, 
such as that of Ref. 0, is necessary to get the correct 
hierarchy in the absence of matter. Moreover, a bulk po¬ 
tential for the radion (the metric degree of freedom as¬ 
sociated with the fifth dimension) is necessary to achieve 
solvable junction conditions when matter is present on 
the boundaries 0. In this case, cosmological solutions 
to order p/K 0 are not sufficient to grasp the particular 
features of the background metric evolution originated 
by the extra dimension and one needs to investigate the 
effect of terms of order {p/VY (as was done in Ref. [Tn|l 
or higher. 

The aim of this article is to go beyond the first order 
approximation in brane cosmology for RS I models with 
two branes. Our approach will differ from Ref. M in 
that we do not consider a bulk scalar field to stabilize 
the radion but include an effective stabili zing potential 
directly into the equations (see also Ref. [ij)- Gonse- 
quently, our perturbative expansion is around the RS I 
solution. The calculations are then carried out in or¬ 
der to show how p^ contributions to the four-dimensional 
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Hubble parameter may affect the model (or may be un¬ 
observable). Such terms are expected as fingerprints of 
the fifth dimension in analogy with the single brane RS II 
framework. The latter case will also be studied as RS I in 
the limit when the distance between the branes diverges. 
Some hints about the possibility of an accelerated expan¬ 
sion driven by exotic fluids with pressure p = w p and 
w > 0 will be presented, thus suggesting the necessity to 
go beyond the second order approximation. 

The paper is organized as follows: in Section m we 
present the complete setup of the model under consider¬ 
ation; in Section ITTTI the second order ansatz is described 
and Einstein equations are perturbatively solved; in Sec- 
tionl^cosmological consequences of the solutions are an¬ 
alyzed and compared to the known brane-world solutions; 
in Section m the analysis of the approximations is per¬ 
formed and, finally, in Section ED some conclusions are 
drawn. For the five-dimensional metric pab we shall use 
the signature (-I-, —, —, —, —), so that g = det( 5 AB) > 0 . 

II. EINSTEIN EQUATIONS 


and /Istab is the Lagrangian of the stabilizing field (for a 
scalar field, see e.g. Ref. 

We shall consider the particular case in which the 
stabilizing mechanism can be effectively described by a 
harmonic potential for the radion, with an effective La¬ 
grangian of the form 

Ces = -V 9 ‘^Hb-bof = -^U{b), ( 5 ) 

when the metric is written as in Eq. o, and the poten¬ 
tial U depends on the component 344. The two 3 -branes 
have opposite tensions and their contribution to the total 
energy-momentum tensor is given by 

rr,A _ b{y- Vi) 

xdiag (Uj -k Pi, Vi - Pi, Vi - p„ V - p„ 0) , (6) 

where i = p, n, and Up = 0 {yn = 1/2) is the position 
of the positive (negative) tension brane. The Einstein 
equations in the bulk. 

Gab = Tab , ( 7 ) 


Let us consider a RS I model perturbed by the presence 
of matter on the two branes. The bulk metric is given by 

ds^ = pAB dx"^ 

= n^{y, t) (V? — a?{y, t) dx* da;* — b^{y, t) dy^ .(1) 


The Einstein tensor for this metric is 



form a system of four differential equations for the three 
independent functions fa = {n,a,b). On using the 
Bianchi identity aG^^ = 0 , it is then straightforward 
to show that the three equations 

( Goo = Tqo 

< Go 4 = 0 (8) 
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where a prime denotes a derivative with respect to y and 
a dot a derivative with respect to the universal time t. 
The energy-momentum tensor in the bulk is that of an 
anti-de Sitter space with the addition of a term gener¬ 
ated by a field which serves the purpose of stabilizing 
the distance between the two branes of the RS I model, 

r 4 = Ag 1 j+fi, ( 3 ) 


are independent and form a complete set. This means 
that a solution to Eqs. © also solves the full set 0 . 
Moreover, since Gab is trivially conserved because of 
the Bianchi identities, the tensor Tab is also automat¬ 
ically conserved regardless of the fact that the effective 
potential in Eq. 0 does not appear covariant. 

For computational purposes, it is convenient to do 
some further manipulation. The first of Eqs. m can be 
replaced by 


F\y,t) 



with 


F{y,t) 


{a a')2 



( 9 ) 


( 10 ) 


which, on intergating along the extra dimension, can be 
written as 


where, as usual, 
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where c{t) is related to the boundary conditions at the 
branes. The conservation equation = 0 yields 

i'= ('!!: +3^1 + 2^) , (12) 

\n a 0 J 


In order to keep track of the various orders in the above 
expansion, it is useful to introduce an expansion param¬ 
eter e by replacing pi ^ epi (and setting e = 1 at the 
end of the computation). We make the following ansatz 
for the metric, 


which is identically satisfied by the solutions of the sys¬ 
tem Instead of the three Eqs. 0 , we shall therefore 
solve the equivalent system 


b'^a^ 


■fa^ (r°o)'d?/ = ^5(t) 


^ b' = {bo-b) (^+3^ + 2^) (13) 

^ G 44 = T 44 . 


Moreover, bulk solutions must satisfy the boundary 
equations given by the junction conditions on the two 
branes. 
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where we assumed an equation of state for the vacuum 
perturbations of the form pi = Wi pi. When pi ^ 0 the 
RS I solution is fully recovered and one finds the usual 
warped static metric with c(t) = 0 and 

nRsiy,t) = aRsiy,t) = exp (-m 60 |2/|) 

(15) 

bRS = bo , 


which also require the well known fine-tuning 


y __y 

Vp — — ^2 ’ 


A = - 


6 mr 

1 ^ 


(16) 


III. THE LOW DENSITY EXPANSION 


n{y, t) = exp {-m bo \y\) [1 -h Sfn{y, t)] (19a) 


a{y,t)=ah{t) exp{-mbo\y\) [l + 6fa{y,t)] (19b) 

biy,t) =bo + 6fb{y,t) , (19c) 

so that the homogeneous scale factor ah{t) is factored 
out, and bo is the equilibrium point corresponding to the 
RS I model [see Eq. (tT^ above]. The solutions to the 
equations can then be completely expressed in terms 
of the functions Sfa and Hh = an I ah, which we expand 
to second order in e as 



We also expand c = c/af as 


c(t) ~ e pp + 

+e^ ( 4 ^^ Pp + 4^^ pI + Pp Pn) ■ (20c) 

Note that the time dependence is just carried by the func¬ 
tions Pi = Pi{t) and that, of all the coefficients appearing 
above, only the fill's in Eq. 12()all depend on y, whereas 
the others are constant. 

In order to proceed with the perturbative expansion, 
one also needs to expand pi. From the conservation equa¬ 
tion 


Pi =-2,H{yi,t) {l + Wi) Pi , (21) 


A perturbative approach can be adopted in brane cos¬ 
mology to investigate solutions to the Einstein equations 
by taking as a starting point the static RS I metric with 
Pi = Q reviewed in the previous section. In fact, in the 
low density regime 




< 1 


(17) 


one can express the corrections to the solution llTHll to 
all orders in pijVi by assuming that the metric functions 
fa = {n, a, b) can be written as 


fa = fRS + Sfa 


(18) 


with Sfa 


'^3 

Pi Pj 


it immediately follows that the time evolution of the mat¬ 
ter densities is adiabatic, since H = a/a ^ and, 
therefore, 




If we now assume that, to second order in e. 


( 22 ) 



(23) 



the coefficients of the above expansion can be related to 
the corresponding ones in Eqs. 120bll and (I2()all for a = o 
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by equating the two expressions for at y = yi up to 
second order, 






S'faiyi,t) 

1 + 6 fa{yi,t) 


2 


With the help of Eq. 121II . one finally obtains 


(24) 


equations ns. This will allow us to determine explic¬ 
itly the coefficients 5fa once the boundary conditions are 
imposed. Let us begin with first order equations. 

At order e, the constraint 113 reads 


f(i)' 


( 1 )', 


/L Pp+fin Pn-^mk 


fW Ai) 

•/6,p rp ^ Jb,n 


= 0 . (26) 


hl^\y) = (25a) 

hf\y) = hfl.-p>{l + w{)h^l;lf^^}{y) (25b) 

-6(1 -b Wn) hill fl%y) ■ (25c) 

A. First order results 

In order to solve the bulk equations order by order, one 
has to substitute the previous expansion in the dynamical 


If we allow pp and pn to be arbitrary functions of the 
time, the above equation splits into two independent 
equations for the coefficients of pi, and one finds 

fb]i (y) = exp (4 mboy) , (27) 

where are constant coefficients to be determined. 

The functions fl^jiyfs can now be calculated by solv¬ 
ing Eq. HU- Since the contribution of the stabilizing po¬ 
tential vanishes at order e, the integral-differential equa¬ 
tion EJ becomes the first order linear differential equa¬ 
tion 


Yl e- 2 p, ( 2 ^ e'^'^^°y-bo cf^ o'*^ -b 5o + 2 m 0 '^ 

i 

On solving for each coefficient of pi independently, one obtains 

r(l) _ ^ „2m bo y „2 mbo y ( fi _ o 1 , „(1) 
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Finally, one can expand the equation G 44 = T 44 to first order. 


(ebohll + 95ow,/ig - 95ocf^e^™^«^ - -b dfcVbgfof= 0 , (30) 


and solve the two equations for (y). The result is 


f(l) _ £ p2mfao y 
J n.i 2 ^ 
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We are now left with six numerical coefficients 


the discontinuity constraints for a'{y,t) at y = yi imply 




b(l) 


„( 1 ) 


(32a) 


and four integration constants 





(32b) 


In order to fix the above, one has to use the junction 
conditions. These four conditions, written in terms of 
the coefficients of pi, form a system of eight equations: 


Pp (3 - 3 h^l -b TO^ -b 3 p„ - hll^j = 0 


^ 3pp(cWe™^o-/bW) 

-bpn ( 3 — fc^TO^ = 0 , 


(33) 
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whose solution is given by 
( 1 ) 

iT'h ^ — 


h,p 3(gmbo_l)’ 

m 


c« = 


p 3(e"i6o_i)’ 


/I \ Ic^ T71 A 

3(emfco-i) 

cW = f . (34b) 

” 3(e™*'«-l) 


The two analogous constraints for n'(jj, t) are both equiv¬ 
alent to the equation 


bo 


m{3wp — 1) + 4uj‘^ bob^^ — l) Pp (35) 


m{3wn-l) + Auj'^bo 6 ^')e™(e™- l) 


Pn — 0 


We see that the junction conditions are not sufficient to 
determine the integration constants 132bll . Such freedom 
is in fact related to the gauge freedom in the choice of 
the initial value for the scale factor and time variable. 
Without loss of generality, and to simplify the second 
order calculations, we then set = 0. The values of 

the c^')’s are related to the choice of the time variable. 
Since one usually considers the negative tension brane 
in RS I as the four-dimensional “visible Universe”, it is 
natural to use the proper time r on this brane and choose 
the c^')’s so as to have n(y„, r) = 1. This can be achieved 
by setting 


which yields 


p 


m(3wp — 1 ) 
iuj^bo (e™^o - 1) 


( 1 ) _ m(3w„ - l)e 

" ~ 4w2 6o (e™*"'- 1) 


(36) 


c(i) = _ 
n,p 


g2mbo 1^3 777,2 ^3^^ — 1) -f 2 6 q (3 Wp + 2)] 

48 m 6 q (e™ — 1) 


c 


( 1 ) 

n,n 


embo [3m^(3wn-l) + 2b^k^uj^ (l-2e-'"'’») (3w„-h 2)] 
48m6ga;^ (e^&o _ i) 


(37) 


and defining the time coordinate r as 

dr = exp ■ (38) 

With this choice, the cosmological Friedmann equations 
can be easily compared to standard ones. 

Let us now comment on the first order results. As far 
as the radion perturbation is concerned, we found 

g4mboy 

<j/b = —— 5 — [(1 - 3wp) Pp + e“™''‘’(l - 3wn)pn] e 
4 Oq 

+0(6^) , (39) 

which was expected, as it is due to the known coupling of 
the radion with the trace of the energy-momentum tensor 
of brane matter. Traceless fluids, such as radiation, have 
no first order effect on the excitation of the radion. If one 
fills the branes with some pressureless fluid, the distance 
between the two branes grows. This effect, being counter¬ 
intuitive for the attractive nature of Newtonian gravity. 


is in fact a consequence of the form of the stabilizing 
potential. Its non trivial contribution to the bulk energy- 
momentum tensor at order e is 

ni^-2bW5h. (40) 

The first order expressions are identical to those for a 
static solution, as they can be obtained by neglecting 
Pi. Every kind of matter on the branes thus acts so as 
to detune the brane tensions from the bulk cosmological 
constant and can be balanced by some constant pressure 
along the fifth dimension. Such pressure is given by the 
first order contribution of (EIHl which increases when dfb 
decreases. 

Note that (TTHll is proportional to the inverse of ui^, 
which represents the effective spring constant coming 
from some stabilization mechanism. When such a con¬ 
stant diverges, the correction Sfb vanishes and the length 
of the fifth dimension is fixed as expected, even if there is 
a finite, tu-independent T 44 pressure term. On the other 
hand, the correction to the scale factor, 
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p 2 m ho y 

- 1) ” - "°'~ °) + P"- - "•'I P> 

+6-"’*" [2tJt^ij*(e^’"‘"''- 2)+3m’(3io„ - l)e*"’'’<"']p„}e + 0(e*) , (41) 

never vanishes when matter is present on the branes, even if that is trace-less. Furthermore a finite, non vanishing 
6 fa can be obtained in the limit of infinite spring constant, regardless of the matter equation of state. 

The carrection to the lapse function, 


p2 m ho y 

= 48m6g02Me-^°-l) ^ '’O _ 32 ho y) + 3 ^2 (3^^ _ 1) p2 ho y ] ( 42 ) 

+e-™'’“ [2blk^u;^3wn + 2) (2 - e^™^"^) + 3m2(3n;„ - 1 )p„} e + ec« Py + ec« + O(e^) , 

I- 


does not vanish when the branes are filled with trace-less 
matter. Note, however, that a vanishing correction can 
be obtained for some exotic fluid with wi = —2/3 and 
negative pressure in the limit a; —> 00 . Apart from these 
exceptions, one has non negligible corrections everywhere 
in the bulk. 

In order to compare the first order results with the 
RS II case of a single brane, we must instead use the 
proper time on the positive tension brane. This is 
achieved by setting 

[3 m?(3wp — 1) -f 2 6 q (1 — 2e™*'°) (3wp -I- 2)] 

dSm^Qo;^ (e"*^“ — 1 ) 

(43) 

[ 3 to^(3w„ — 1)-I-25q (3u>„-I-2)] 

48m6gti;^ (1 — e™^«) ’ 

and letting 60 ^ 00 . The Friedmann equation is simply 
obtained by keeping 0{e) terms in Since the 

first order four-dimensional Hubble parameter is homo¬ 
geneous, it reacts to all the sources along the y direction. 
On the positive tension brane one has, to first order in e, 

2 TOfc 2 (e™'’Vp + e”™^Vn) ho^oomfc^ 

- z -71-e —^ — 1 —Ppei(44) 

P 3 (em 6 o_x) 3 rp / 

regardless of the value of w. When the negative tension 
brane is moved to infinity, its contribution goes to zero 


and one recovers the usual first order effect in brane cos¬ 
mology. 

On the visible brane, the Friedmann equation is 
slightly modified by the rescaled time parameter r, 

2 ^ wfc^ (e^™'’Vp + Pn) Pp^o mP pne 
n- 3 (e"4 bo _ 1 ) ' 3 (e"! bo _ 1 ) ■ 1 '' 

The acceleration equation is homogeneous as the Fried¬ 
mann equation and becomes 

a{yi,t) _mk‘^ [e""'’'>(l -f Wp)pp + e“™'’“(I -f Wn)pn] 
a{yi,f) 6 (l-e™*'«) 

(46) 

which has the weighted brane fluid energy densities as 
sources. 

B. Second order results 

We are now ready to evaluate 0{e^) corrections to the 
vacuum solution RS I. The procedure will be analogous to 
the one used for first order results in the previous section. 

As a first step, one can impose the constraint 03 in 
order to find the dependence on y of the second order 
coefficients in Sft- We are then left with three inhomo¬ 
geneous equations obtained by setting to zero the coeffi¬ 
cients of the independent matter densities in 


,(2)' 2 , f{2)' 2 , --(2)' . , f f{2) 2 , j(2) 2 , f{2) \ 

fb,p Pp + fb,n Pn + flrn PpPn-imbo / p^ -f /// -P Pp j 
m (1 - 3 Wp)2e®n 2 , 

24 5ga;4 (e'^^o _ l) T 
m(l-3w„)2e2™*'« 


" (6 + bl k^ - bl k^ '>»] 

24 5ga;4 (e'^^o - 1 ) ^ o j o j 


Pp 

2 


24 5gw4 {e^bo _ 1) 
m(l- 3u;p)(l - 3w„)e'"'’o 
24 6ga;4 (e™'>o - 1 ) 


[e^ (6 m^ + bl bl k^ bj^] pI 

|‘2e2mboy ^0^2 _|_ ^2 ^2 _ ^2 ^2 ^2 ^4 _|_ g^bo^j — Q 


(47) 
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which contains the hrst order parameters previously determined. The solutions are 

(1 - 


/S = - 


96 6qu; 4 (e™^o — 1) 


6 TO^ + 6 ^ - 2 hi e™ 


(48a) 

(48b) 


/ON f , . _ ‘? 7 /t p6m6o(y-l) ) 

Jtl=‘'"‘^'' hg) - (6„. + tg (I +e'"<.)] I ,(48c) 

( 48o5w‘^ (e™''« — 1) J 


where 


6 ( 2 ) fo( 2 ) 5 ( 2 ) 


(49) 


and nine parameters related to the radion, the four¬ 
dimensional Hubble parameter and c(t) respectively, 


are integration constants to be determined from the junc¬ 
tion conditions. Once we plug the into the second 

order terms in Eq. m, we get the equations for the 
( 2 ,) 

which are not displayed for the sake of brevity. 
Finally, by solving G 44 = T 44 one obtains the cor- 

rections /„ (’s. The results will contain six integration 
constants from the solutions of the first order differential 
equations for and 


c 


( 2 ) 

a ,2 ’ 



(50) 


h 


( 2 ) 
h,i ’ 




(51) 


where i runs over p, n and m for second order quantities. 
Analogously to the first order case, one can fix the coef¬ 
ficients eu by imposing the junction conditions, which 
form a system of nine independent equations. Nonethe¬ 
less, one is again not able to fix the constants CT . Here 
are the solutions obtained at the end of the calculations 
described above: 


= -2 \3m^{l - 3wp f-8 blk^oj^] (6m^ - blk^uj^) 

4608m2 64a;4 (e™"-"-!) 0 JV 0 J 

-f8 6^/c2w2e'"*'«(®*'+i( [11 to 2 (3wp - 1)-k 2 52 fc2a;2] 

-4 bl k^ uj^ e™ (2?/+i) [3 ^2 (-g ^2 ig _ 7) g 52 ^2^2j 

-f4 6^A:2^2g2mbo(y-ri) [27 to2 + 2wp - 1) + 16blk^uj^] + 

_l_g8mhoy 1^45_ 3i(;p)2 _|_ Zrn‘^bQk^uj‘^{9Wp — A2wp + 13) — 46gA:"‘tt’"^] 

_2g2mho (2y+i) [gi^4(9^3 g^2 _ {27+ 42wp - 1) -f 86^A:^a;^] 

— 2e'"^“ (4y-|-l) j-g^ ^4^g^3 _|_ g^2 _ _ 3 774,2^2^2^2^9^2 _|_ 0^^ _|_ 29^ _|_ 45 &4^4^4j | _|_ ^^2^ (52a) 

= -2 (2e2™^'>(2'+i( r3m2(l - 3u;„)2 - 86gfc2u;21 (6m2 - blk'^Lo^) 

4608m254,,4 (embo _ 1)2 I L ^ o J V o J 

[llm2(34(;„ — 1) + 26ofc^w2] — 46gA:2a;2e2™*'“ [3m^{9w^ + 18w„ — 7) -I- 86Qfc2a;2] 

-k46ofc2w2e™^“ [27 to2(3u;^ -k 2wn - 1) -k 166ofc2u;2] 

_|_g6mboy 1^457774^4 _ JiWnf + 3w?bQk^uj^ {9vj^ — 42u>„ -|- 13) — 460^^0;'^] 

_2e2mboy [gi7774(97g3 g^2 _ _ 4) 3m%lk^oj^{27wl + 42w„ - 1) -f 86ofc‘‘w'‘] 

_2gmbo(2y+i) [347774^9^3 9^2 _ 77^^ _ 4) _ 2,m%lk^oj'^{9wl + 6w„ -k 29) -f m^k'^Lo'^] | -k 0^2^ (52b) 
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„m 6 o( 2 y-l) 

f^}n = -^- 2 ~ — Bwn + 11 ) + 26QA:^a;^l 

2304m2&X (e’”'’“ - 1) L ^ p A p ; 0 J 

+46gfc^a;^e'"*'“^^*'''"^^ [m^(3'u;„ — l)( 6 w„ — Gwp + 11 ) + 26gfc^w^] 

+46gfc^a;^e'"*'“ [3m^(18Wp — 9wpWn + 9wp — 9w„ — 1 ) + 4&Qfc^a;^] 

+46gfc^a;^e^’"^° [3m^(18'u;^ — 9wpWn + 9wn — 9wp — 1 ) + 46gfc^w^] 

_l_g 6 m 6 oy \^ii^rn^{ 2 ,yo^ — l)(3t(;„ — 1) + 3m‘^bQk^uj^{9wpWn — 21wp — 21w„ + 13) — 4&gfc‘*a;"‘] 

_g 2 m 6 oy 1 ^ 977 ^ 4 ^ 3 ^^ _ _|_ Qyj^yj^ -(- 27Wp + 3w„ + 11 ) 

+‘iQm'^b‘^k^uj'^{iw^ — 3wpWn +Wp — 2wn — 1 ) + 86 gfc‘*w‘*] 

_g2m6o(y+l) [ 9777 , 4 ^ 3 ^^ _ l)(lg 77;2 _|_ Qyj^yj^ -\- 27Wn + 3w„ + 11) 

+36m^bQk^uj^(3w^ — 3wpWn +Wn — 2wp — 1 ) + 86 gA:^a;^] 

_gmho(2y+l) |^ 8 l 7774^0^3 _|_ 0^3 _|_ _|_ 2,WpW^ + 7Wp + 7w^ + ^WpWn — Wp — Wn — 2) 

+6m^&Qfc^w^(18Wp + 18w^ — 9wpWn + 21wp + 2lwn — 1) + 166gfc^w'^] } + (52c) 


p2m6oy 

f^l = -2 |—86ofc^w^(3wp + 2)e^™ ° r3m^(3'u;p — 1) + 2b1k'^uj‘^{3wp + 2)1 

4608m264^4 (emho _ 1)2 I o v p ; L v p ^ o v p 

-86gA:2w2em6o(4y+i) [7772(31^3 130^2 2^Wp - 31) + 2blk‘^uP{<owp + 5)] 

j^2e'^mboy ^0^2 _ ^2^2^2^ [Zm^{'iwp - 1)^ - ^blk^uj"^{9wl + 6wp - 1)] 

—46gfc^a;^e"‘*'“ [—3m‘^{63Wp + 30wp — 17) + 46gfc^a;^(9Wp + 6wp — 1)] 

+46gA:^a;^ [—27m^(18Wp + 27Wp + 4wp — 5) + 86gfc^w^(9Wp + 6wp — 1)] 

_|_g6m6oy 1^457774('377;^ _ 1)2 _|_ 37772^2^2^2)111.^2 _|_ 3077;^ _ 23) -|- 46Qfc‘*W^(6'U;p + 5)] 

_2gm6o(2y+i) [sim^{9wl + 9^^ _ 77;^ _ 1) _ SmHlk"^iv"^{162wl + 279wl + 60wp - 25) 
+8b^k‘^uj^i9wl + 6wp - 1)] + [-9771^(81^3 + 90^^ - 15wp - 8) 

+3m^6oA:^w^(1627(;3 + 2437Pp + I2wp — 53) + 46gfc‘*w‘*(9Wp + 2417;^ + 14)] } + (52d) 


p 2 mho(y-l) 

" 4608m2bX(e^5o-l)^ {-45gfc2u;2(67c„ + 5)6™*-° [27m\3wl + 2wr. - 1 ) + Sb^k^u;^] 

— 86 oA:^w^e^'"*'°^ [to^(81w 3 + 180w^ + 24w„ — 31) + 2blk^uj‘^{6wn + 5)] 

+460^:^07^6^'"*'” [3to^(45w^ + 2Awn — 13) + 45Qfc^o;^(6w„ + 5)] 

_ 2 embo( 2 y+i) [ 8 l 7774 ( 9 w 3 + Qw^ _ 77 ;^ _ 1 ) _ Sm^blk^u^{162wl + 243wl + 48w„ - 17) 
-8blk^uj^{Qwn + 5)] + 2e2'"'’°(*'+i) [9r77^(3w„ - 1 )^ 

-3m%lk^uj^{^3wl + 30w„ + 1) - 45;]fcV(6w„ + 5)] + e®'"*’”*' [45r77'‘(3w„ - 1)^ 
+3 to^6qA:^07^(117w^ + 30w„ — 23) + 46 qA:^07^(6w„ + 5)] + 2e^'"*'”^ [—81 to^(9w3 + 9w^ — «;„ — 1) 
+3to^6qA:^07^(162w 3 + 243w^ + 12w„ - 53) + 86gfc'‘o;'‘(6w„ + 5)] } + 


(52e) 
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f(2) 


-^ \Ablk^uj^{3wp + [3771^(1 - 3w„) + 2blk^uj^(3wn + 2)1 

2304m264^4 (gmbo _ 1)2 I 0 ' P ' L ^ 7 0V 7J 

—4:blk^uj^e^"^’’°^ [m^{3wp — l)(18Wp + %WpWn + 45wp + 24w„ + 31) + 2blk^u;'^{3wp + 3wn + 5)] 

—4&Qfc^w^e’"^‘’^®^''-i) ^rn?{3wn — l){18w‘^ + 9'WpW„ + 45'u;„ + 24wp + 31) + 2blk‘^uj^{3wp + 3wn + 5)] 

—1^377^2^3^;^ _|_ 2)(18io^ + 9wpWn — 12wp + 27w„ — 8) 

+2&Qfc^a;^(18'u;p7(;„ + ISwp + 15w„ + 13)] — 46g/c^a;^e™*'“^^^''-i) ^—2blk‘^Lo‘^(9wpWn + 3wp + 3wn — 1) 
+3m^(37i;p + 2)(18Wp + 9wpWn + 27wp — 15w„ — 7)] + [45m‘*(3wp — l)(3w„ — 1) 

+3m^6Qfc^a;^(117u;pW„ + 15wp + 15w„ — 23) + 4&QA:^a;^(3wp + 3wn + 5)] 

_|_gm6o(4y+l) |]_gl77j4^0^3 _|_ 0^3 _|_ _|_ ZWpW^ + 7w^ + 7w^ + AWpWn — Wp — Wn — 2) 

+6m^blk'^u}'^{54Wp + Mw\ + 27w^Wn + 27wpw'^ + 117Wp + 117w^ - 9wpWn - Qwp + 12u;„ - 49) 
+8&Qfc4a;4^g^^y^^ _l_ i2iUp + 12wn + 14)] + [—9m4(3wp — l)(18Wp + 9wpWn + 27wp + 3wn + 11) 

+54m^6gfc^a;^(6Wp + SWpWn + 13Wp — 2wpWn + ^Wp — 3wn — 1) — 4&Qfc‘*w‘*(9wpW„ + 3wp + 3w„ — 1)] 

+e‘^mbo(2v+i) 2wn + l)(34i;2 + 2zi;„ - 1) 

+6m^6gfc^w^(54w^ + 27wpw'^ + I17wn — Ihwp + 30w„ — 13) 

-AbQk‘^Lo'^{9WpWn + 3Wp + 3Wn - 1)] } + C^n}m ■ (521) 


IV. SECOND ORDER COSMOLOGY 


A. RS I 


Let us now look at the physical consequences of the 
above results as applied to two different cosmological sce¬ 
narios. The Friedmann and acceleration equations will 
be showed and their phenomenology investigated in dif¬ 
ferent regimes of to and for different matter equations of 
state on the branes. 


We begin from the case with two branes HI for which 
we shall study how cosmology would be described by ob¬ 
servers on the negative brane. In order to achieve that, 
one can use the corresponding proper time by fixing the 
as in Eq. e3, choosing 


g 2 m 6 o 

4^], = - 5 - {9m'^{3wp — 1) \3e'^"^'^°{3wl + 4wp -I-1) -I- 18e'"^“(3w^ -I- Awp -f 1) — 12wp + 4] 

4608m2&4,,4 (gmho _ 1)2 f ^ ^ L ^ P P ' '' P P ’ P J 

-ni^blk^uJ^ [e^^^°{32Awl + 153wl - 102?«p - 91) - 6e'"^«(162u>3 -t 207?w^ -k 12u>p - 65) 

-k54(5u>p -k 2wp - 3)] -t 46;)fcV [3e^”^’’°{6wl -t 6 wp -h 1) -t 2 (l - 2e™'’°) (9^2 -t 6wp - 1)] } (53a) 

= 4608m264^4\en.feo _ i )2 " 1 ) - 1 ) + 18 (e-^« + l) {3wl + 4n;„ + 1 )] 

- 771 ^ 62 ^ 2^2 [e2'"''0(81w2 _ 90w„ - 75) -k 2e'^'’°{162wl + 279wl + 192w„ - 5) 

-6(162w 3 -k 2A3wl + 607(;„ - 37)] -t Ab^k^u;'^ (e™^“ - 2^ ( 6 w„ -k 5)| (53b) 

pmho 

c(i]L = , , —7 -^ {97774 [e2'"^“(3?7;„-l)(18w2+ 977 ; u;„_6?Cp-k 27w„-t 14) 

2304m26;)a;4 (e'"^o - 1) 

9e™'>0 _|_ 0^3 _|_ 2,w^Wn + 3wpw‘^ + 7Wp + 7?«^ -I- 4wpWn - Wp - Wn - 2) + {3wp - l)(187(;p 

-|-9wp7c„ -I- 27wp + 3wn + 11)] — m^b^k^uj'^ (lOSw^ -I- 547Cp7c^ -|- 2347c^ — 9wpWn + 87wp 
—39wn -I- 1) -I- 2e™*'°(54uip — 16277;^ -|- 27WpWn — 81WpWn + 1177Cp — 35177;^ — ibWpWn + 30wp 
— 96Wn + 11) — 6(54771^ -I- 27WpWn + 11777;^ — 36WpWn + 30Wp — 33Wn — 19)] 

—AbQk'^uj^ [(e^™^° — 2e'"*'“) (377;p -I- 2){3wn -|- 2) -|- 9wpW„ + 3wp + 3wn — l] } 


(53c) 
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and finally rescaling n{y,t) to satisfy the condition n(l/2,T) = 1. The Friedmann equation we are interested in is 
given by the second order expression of the Hubble parameter at j/ = 1/2 as a function of r, 


H2(i/2,r) = 




3 (^mbo _ 1) 


{pn + - iwn + 1) (3Wn - 1) 


{0^^° + 1) 
3262u;2 {ernbo _ 1) 


2 Pn 


“hi') 1 

— + [l8wp - 9wp{wn - 1) - 9wn “ 1 + 2e™''“(9w^ + 9u;„ - 4)] ^ e™'’“ pp p„ 


486gw2 (e’"''') - 1) 

' m2fc2(3wp-1) (e™'’« + 1) 


36 


966ga;2 (e'"*'o - 1) 


[3u;p + 7 - 4e"^'''>(3«;p + 4)] I e^™^" . 


(54) 


Note that the results (I53all - (I54II do not depend on the integration constants c^^]’s and c^^]’s, which reflects the fact 
that the three-dimensional spatial curvature has been set to zero ab initio. 

The Friedmann equation contains coefficients up to second order in the vacuum perturbations of both branes. First 
and second order contributions in Eq. (IHUi are consequences of the adiabatic regime of the five-dimensional dynamics 
which determines the value of the integration constant c(t) in Eq. (Hill . The value of the latter is affected by the 
presence of matter on both branes through the junction conditions and, for instance, up to 0 {e^) is given by 


lim 

UJ —»-oo 


lim 

CJ—^■OO 


mP {pp + e”^^°pn) 
3 (e'"*'o - 1) 


fc4 -p2e™^« - 1 n /e'"'’“-fl\^ -h 2e-™'’« - 1 , 


(55a) 

(55b) 


in the limit of infinite spring constant. Furthermore, as 
previously noted, this effect is also a consequence of the 
radion field potential acting as a source in Eq. One 

can expect to cancel some terms in Eq. (ISH) by arbitrar¬ 
ily increasing the spring constant of the effective radion 
potential in order to decrease the radion shift from equi¬ 
librium. This mechanism partially works as the first or¬ 
der contributions to the bulk potential vanish along the 
time-time direction, whereas only second order terms, 
which depend on the matter equation of state, cancel 
when w ^ oo. One is then left with 


lim H^{1/2,t) 


mk'^ (pn + 


3 (e”*'’o - 1) 


_ p-mbo 

36 


fc4 
36 ' 


pI :(56) 


which is analogous to what has been obtained in Ref. M- 
The matter on the positive tension brane appears at sec¬ 
ond order with the role of some “dark” fluid and acts as 
a sort of Brans-Dicke field which adiabatically modifies 
the Newton constant perceived in the visible Universe. 
This behavior is somehow inherited by the dynamics of 


the radion which is known to modulate the strength of 
gravity on the visible brane. 

Note that while first order coefficients are positive def¬ 
inite, irrespective of the brane tension, the sign of sec¬ 
ond order ones depends on the matter equation of state. 
Differently enough from unstabilized brane cosmology, 
this fact implies that leading order cosmological equa¬ 
tions have the correct behavior on both branes. Letting 
the energy density pp ^ 0 in Eq. (IMIl . one obtains 


772(1/2, r) 


{Wn + 1) (1 


3Wn) 


3m^ fc2 (e"^bo _j_ 1^ ^ 

32 6§w2 (e'"*'o - 1)^ 


w fc2 p„ 

+ 3 (embo _ 1) > 


(57) 


which has the usual first order solution for both radiation 
and a cosmological constant on the visible brane. A mat¬ 
ter dominated Universe would otherwise generate second 
order corrections. 

We now come to the equation for the acceleration, 
which has the general form 
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a(l/2,r) 


z(l/2, r) 6(1 —e’"*'o) 


[(3w„ + l)pn + + l)/Op] + {Wn + 1) (3w„ - 1) (3w„ + 2) 


3m^k'^ (e^^o + 1^ 
3262u;2 (e^feo _ 1) 


2 r n 


~ i 144^^^" + 1) + 


966oW^ (e™^o _ 1)" 


[e^^^°{3wn — l)(18w^ + 9wpWn + 15wp + 39w„ + 23) 


+9e™*'“(6Wp + Wn + 3WpWn + Swpwl^ + llWp + llw^ + 2wpWn - 4) 
-\-Mwp — 27wpWn + 99wp — ISwpWn + 2Awp — 13] | e™*'“ pp pn 
2 m?k‘^{Zwp — 1) + l) 


k^ 

+ <---(3i(;p + 1)^ + 2 

1 144 966qW^ (e^&o _ i) 


[e'"^“(27w2 + 54u>p + 23) - 9wp - 5] ) e^’^^'Vp , (58) 


and, for cu ^ oo, reduces to 
d(l/2,r) 


m 


i(l/2, r) 6(1 — 


[(3Wn + l)pn + e^ "* {3Wp + 1) Pp] 


~ 144^^^” + l)(3Wp + 1)0™^“ PpPn + 1)^6^'"*'“ Pp , 


(59) 


which is again analogous to the result of Ref. M- The coefficient of Pp in Eq. is positive or zero in this limit and 
provides an accelerating contribution to the equation. The coefficient of the mixed term has a positive value when 
just one fluid has Wi < —1/3. For pp ^ 0, one is further left with 


d(l/2,T) mk^{3wn + l) 
a(l/2,T) 6(1 —e™'’o) 


H~ 1 ') 

Pn + {Wn + 1) (3ZC„ — 1) (iWn + 2) — - —^ p„ , 

326§a;2 (e^^o - 1) 


(60) 


where the second order contribution is inversely propor¬ 
tional to Lo and vanishes for radiation and a cosmological 
constant. This peculiarity leads to the standard cosmo¬ 
logical evolution up to 0 {p^) until the matter dominated 
era. 

A singularity in the lapse function n{y, t) for bo ^ oo 
prevents us from analyzing the correct limit when the 
distance between the branes becomes infinite, we shall 
thus comment on this problem in the next subsection. 
On setting pp 0, Eq. 15611 admits the finite but trivial 
limit 


This result is due to the reduced strength of the gravita¬ 
tional interaction at infinity. 


B. RS II 


lim H^{1/2,t) = 0 , 

bo — >-(x> 


in which one also has 


lim 

bo—*oo 


a(l/2,T) 

a(l/2,r) 


= 0 . 


' ^ On can think of the RS II model as the limit of RS I 

in which the distance between the two branes becomes 
infinite, thus one expects that only pp contributes in this 
limit. The cosmological proper time is now the one on 
(61b) the Planck brane and is recovered upon choosing 


= --- 2 {9w"‘(3wp - 1) [l8e^™'’“(3wp-I-4?iip-I-1)-I-18e'"*'°(3wp-I-4wp-I-1) 

4608m264^4 ( 3^60 _ i )2 1 ^ p ^ L \ v p / v p 

—21wp -I- 7] -I- m^bQk^uj^ [ 6 e^'"*'“ (162Wp -|- 2A3Wp + 60wp — 37) 

- 2 e™'’° (162^3 -p 279 u> 2 -h 192wp - 5) - 81 -h 90wp -f 75] 

+AbQk‘^uj‘^ (1 — 2e™^°)^ (6wp -f 5)| 


(62a) 
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c(2) = 


-2mbo 


4:608m^bQLL>^ _ i)' 


{9m^(3w„ - 1) [4e2'"^«(l - 3w„) + 18e^^°{3wl + 4u>„ + 1) 


+54^2 + 68wn + 21] - m^blk^uj'^ [Me^"^^°{5wl + 2w„ - 3) 

(162^3 + 207^2 + 12w„ - 65) + 324^3 + lh8wl - 102w„ - 9l] 
+^blk^uj^ [2 (6^™'’° - 2e'"^«) {9wl + 6wn - 1) + 3{6wl + 6u>„ + 1)] } 


(62b) 


c(2) = 

^n.m 


^mbo 


2304to 2&4^4 (gmho _ 1)^ 


{9m'^ [e^'”^“(3w„ — 1)(18?«^ + 9wpWn + 8wp + 27w„ + 11) 


+9e""'’“(6Wp + 6w^ + SWpWn + iWpW^ + 7Wp + 7w1_ + AWpWn - Wp-Wn- 2) 

{3wp — l)(18Wp + 9wpWn + 27wp — 6wn + 14)] + rn^blk'^uj^ [6e^'"*'“(54w^ + +27wpw'^ 

+117ui^ — 36wpWn — 33wp + 30wn — 19) + 2e™*'“(162wp — 54^^ + 81WpWn — 27wpw'^ 

+351uip — 117ui^ + 96wp — 30wn — 11) — 108u>p — bAwpWn — 234:Wp + 9wpWn + 39wp — 87w„ — l] 


—Ab^k^uj^ {9wpWn + 3wp + 3wn — 1) + (l — 2e™^“) (3u>p + 2)(3w„ + 2)] } 
together with Eq. 63- The Friedmann equation, for a hnite bo, reduces to 


(62c) 


H^(0,t) = 




3(gmho _ 1) 


326^^2 ( e ™'>0 - 1)2 


J m^k^ ^Q^nbo _|_ e"^^°{18w'^ — 9wpWn — 9wp + 9wn — l)(9uip + 9wp — 4) 

I 36 246^a;2 (e'"''<i - 1)^ 

i k* m^fc^(e-^° + l)(3m„-l)[e-^»(3u;„ + 7)-4(3u;„ + 4)] ] ^ 

\ 36 966Qa;2 (e™*'o — 1)^ J 


^ — m bo 


Pp Pn 


(63) 


r 


which is similar to Eq. itSl) but has the hnite, non trivial 
limit 


lim H (0,t) = — 

bo—^oo o 


■ Pp 


(64) 


This is precisely the standard Friedmann equation one 
has in four-dimensional cosmology. Furthermore, the 
equation for the acceleration to second order in terms 
of the time on the positive tension brane is given by 


d(0,t) mk'^e^^° 


a{0,t) 6(1 —e'"'’o) 

f k* 


[(3r(;p -I- l)pp -I- {3wn - l)e ^ -I- {wp + 1) (3uip - 1) {3wp + 2) 


(e™''‘> -h l) 
326qu; 2 (e™^o _ i) 


_ 2 

2 Pp 


m 


■fc^ 


1 .. (3^<^p + l)(3u'n + 1)-p 

144 966Qa;2 (e™^o — 1) 


[e^™^“(54r(;^ _l_ 27wpw‘^ + 99w‘^ - 18wpWr, 


-21wp + 2Awn - 13) -I- 9e™'’“(6Wp -I- 6w^ -I- 3w^Wn + 3wpw‘^ + llWp -I- llw^ -I- 2wpWn 
+Wp -I- iCn — 4) -I- (3wp — l)(18wp -I- 9wpWn + 39wp -I- 15?n„ -I- 23)] | e“’"^“ppp„ 


r 


(65) 


and, compatibly with Eq. (IMIl . yields 


a(0, t) 

bo^oc a(0, t) 


lim 


m fc2 


(3iCp -I- 1) Pp 


in the limit of inhnite brane distance. Thus second or¬ 
der effects, typical of brane cosmology, become more and 
more negligible when the distance between the branes 


6 


( 66 ) 



















grows, which consequently leads to unobservable devia¬ 
tions from standard four-dimensional General Relativity. 


V. APPROXIMATION ANALYSIS 


The results obtained so far hold with the assumption 
that Pi M'^ where M is, in general, the natural mass 
scale of the model. By taking all the mass parameters to 
such a natural scale, one expects the solutions to second 
order provide a good approximation for Eq. m- In this 
regime, however, second order effects are certainly sub¬ 
leading and thus insufficient to signihcantly alter hrst 
order behavior. We shall hence present below a numeri¬ 
cal analysis of the validity of our approximate solutions 
in the attempt to extend the range of the parameters in 
which ours results hold valid and widen the conclusions 
one can draw from second order expressions. In particu¬ 
lar, we are interested in possible deviations from standard 
cosmological equations due to terms of O(e^) in the first 
and third of Eqs. ©. Note that all numerical results will 
be obtained by setting the expansion parameter e = 1 as 
previously prescribed. 

In order to test our approximations, since Eqs. are 
not analytically solvable, we substitute the second order 
solutions into the exact Einstein equations to obtain an 
estimate of their non vanishing remainders which we then 
compare with the leading contributions (satisfying the 
corresponding approximate equations). Since the time 
dependence is contained in pi in our expansions (EHall - 
(I2()cll . it is also convenient to trade the time for p. We 
thus divide Eq. m into the following five terms 


Vi{y,P) = 
V2{y,p) = 
V3{y,p) = 
V 4 {y,p) = 
ii5{y,p) = 


n 


_ a'^{y,t) 

‘^{y,t) a?{y,t) 

1 ( a'{y,t) 

b‘^{y,t) \a{y,t) 


2 


6 a^{y,t) 



cjt) 

a^{y,t) ’ 


(67a) 

(67b) 

(67c) 

(67d) 

(67e) 


in which t = t{p) in the right hand sides is understood 
as the time when p = pn (pp will be chosen either equal 
to Pn or zero). The sum. 


Rr,{y,p) , ( 68 ) 

i 


evaluated on the second order solutions yields Rrf{y, p) = 
0{e^) as a measure of the corresponding error. Similarly, 
the third of Eqs. CSl) may be written as the sum of the 
following six terms 
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^i(.y,p) = 
6(2/, p) = 
6 ( 2 /, p) = 

6 ( 2 /, p) = 
6 ( 2 /, p) = 
6 ( 2 /, p) = 


1 a{y,t) 

a{y,t) 

_^ [ a'{y,t) y 

b‘^{y,t) \a{y,t)) 

1 a'{y,t) n'{y,t) 

b'^{y,t) a{y,t) n{y,t) 

1 / a'^iy,t) _ a{y,t)h{y,t) 

\a^{y,t) a{y,t)n{y,t) 

-y(A + £/) 

uj'^ b{y,t) [b{y,t)-bo] , 


(69a) 

(69b) 

(69c) 

^(69d) 

(69e) 

(69f) 


and 


R^{y,p) = Y,^i , (70) 

i 


with R^{y,p) = 0{e^) for the same approximate solu¬ 
tions. One may now assume that approximate metric 
functions computed to O(e^) are accurate approxima¬ 
tions of exact solutions to Eqs. l|I21) if 


| 77 i|>|Ri| and |^/| » |i? 2 | , (71) 

for every term in Eqs. (I67all - 169111 evaluated to O(e^). (It 
should actually be sufficient to satisfy the above condi¬ 
tions for the leading terms of each equation.) Thro ugh- 
out this section, where unspecihed, a natural choice 
of dimensionful parameters is considered. 

Fig-dshows the functions 77 / and evaluated to second 
order in e and the corresponding Rp and R^ for pi = 
2 • 10“^ and Wi = 0.5. The box in the first row on 
the left shows the absolute values of the leading terms 772 
and 773 , which are roughly one order of magnitude larger 
than the remaining terms ry ’s displayed in the plot on the 
right along with Rp (the dotted line). In the second row 
of Fig.[Tl the modula of the leading terms ^ 2 , 6 ,^5 are 

plotted on the left while the remaining coefficients and 
R^{y) are presented on the right. Thanks to the relatively 
large amplitudes of the leading terms, one can rely on the 
second order approximation even when the pi’s are not 
too small, implying that a great improvement over first 
order results can be achieved in a regime where pi ~ 
10“^. In this case the errors produced by truncating the 
expansion to 0 {e^) are much smaller than one percent. 

Let us further consider the particular case pp = 0. This 
choice is made in order to study the cosmological conse¬ 
quences experienced on the negative tension brane (at 
y = 1 / 2 ) generated by second order terms proportional 
to Pn- We therefore use the proper time t on the nega¬ 
tive tension brane M In each graph in Fig. |21 we plot 
the squared Hubble parameter on the negative tension 
brane of Eq. 0 and the acceleration of Eq. lUinil for a 
given energy density pn as a function of its equation of 
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y 



0 0.1 0.2 0.3 0.4 0.5 


y 



0 0.1 0.2 0.3 0.4 0.5 


y 



y 


FIG. 1: The graphs on the left contain the plot of the absolute values of the leading terms among rji (above) and (below) to 
The graphs on the right show the subleading terms among t]i (above) and (below) and the corresponding remainders 
Rri and to C>(e^) (dotted lines). All plots are for pi = 2 ■ 10~^ and Wi = 0.5 and cover all the bulk between the two 
branes. 


state Wn- The plot on the left in the first row is for the 
small density = 10 “^ ^ and shows a behavior 

which is typical of standard four-dimensional cosmology. 
In fact, only Wn < —1/3 leads to an accelerating phase. 
This trend is modified by higher densities, as it emerges 
in the remaining graphs of Fig. [3 In particular when 
Pn = 2 ■ 10 “^ (plot on the left in the second row) the 
second order corrections seem to provide an accelerated 
regime for 0 < Wn < 1. The acceleration appears ampli¬ 
fied when Pn = & ■ 10 “^ (plot on the right in 

the second row) or higher. Note however that the Hub¬ 
ble parameter, a positive definite quantity, constrains the 
region swept by Wn which is not allowed to reach unity. 
Finally the intermediate case is showed in the plot on 
the right the first row: due to the second order effect, 
exhibits a dependence on Wn otherwise not present. 

The four regimes described above have to be tested 
with particular accuracy because the acceleration and 
Hubble parameter plotted in Fig. |21 are not the leading 
terms in Eqs. m and could thus be comparable with 
the remainders. In this case, it is somehow possible that 
the remainders significantly modify the behavior. We 
first note that ~ rji and d/a ~ ^i, as defined above 
in Eqs. (ItiTall and (ItibaH . hence we can use rji and in 
place of and d/a respectively. We then plot in Fig.O 


the ratio between the remainder i?,, (evaluated to second 
order) and the squared Hubble parameter 771 evaluated 
to first and second orders. The four plots show this ratio 
for different choices of and Wn (the same as in Fig. |5| 
and in the same order) as a function of y. In partic¬ 
ular, we choose Wn = 0.95 in order to explore regions 
where the acceleration has an unconventional behavior 
for Pn = 10-2 M4, 10-1 M-i, 2 • 10-1 and Wn = 0.65 
for Pn = & ■ IQ-i Afi. Apart from the last case, the cor¬ 
rections given by the neglected terms cannot significantly 
modify on the negative tension brane. In the first line 
of Fig.^ the ratios between (to second order) and the 

acceleration to first and second order are analogously 
plotted. For low enough energy densities, the second or¬ 
der expressions appear to be good approximations. On 
the other hand, in the second row it is shown that, in 
the unconventional regime of the acceleration (that is 
Pn = 2 ■ IQ-i Afi and Wn = 0.95, and pn = 6 ■ IQ-i 
and Wn = 0.65) the two terms are of the same order of 
magnitude. This shows that one should go beyond 0{e'^) 
in order to determine the true behavior of d/a for such 
equations of state. 



















15 





° 1 -1 - 0.5 0 0.5 1 

Wn Wn 


FIG. 2: Plots of H^{1/2,t) (solid line) and ci(l/2, r)/a(l/2, r) (dotted line) to 0{e^) at a given time, as functions of Wn, for 
pn = 10“^ 10“^ 2 • 10“^ and 6 • 10“^ (from top left to bottom right). 


VI. CONCLUSIONS 

We have computed approximate cosmological solu¬ 
tions of five-dimensional Einstein equations for Randall- 
Sundrum models in the presence of a radion effective po¬ 
tential. The calculations were performed up to the sec¬ 
ond order in the energy densities of the matter on the 
branes and assuming an adiabatic evolution of the sys¬ 
tem. Our approach differs from Ref. M in that we do 
not include a specific bulk field to achieve stabilization, 
and is therefore more general. Interestingly, their results 
are recovered in the limit of very large warp factor and 
radion mass. For the RS I model with matter localized 
only on the negative tension brane, we found negligible 


corrections for the Hubble parameter in the case of ra¬ 
diation or cosmological constant, thus supporting one of 
the main results of Ref. M- For RS II, we found negli¬ 
gible corrections for the equations of state just described 
and in the limit when the distance between the branes is 
taken to infinity. 

On inspecting our results, we finally found some evi¬ 
dence of an accelerating phase for a wider range of val¬ 
ues of the equation of state Pn = Wn Pn on the negative 
tension brane if the distance between the branes is fi¬ 
nite. However, one should then carry the computation 
to higher orders, since such an effect appears near the 
limit of validity of our perturbative expansion. 
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FIG. 3: The two graphs in the first row show the ratio between 7?^ to 0{e^) and rji ~ to 0(e) (dotted line) and to O(e^) 
(solid line) at a given time, for pn = 10“^ and Wn = 0.95 (left) and for p„ = 10“^ and = 0.95 (right). In the two 
graphs in the second row, the same ratios are are given for = 2 • 10“^ and w„ = 0.95 (left) and for pn = 6 ■ 10”^ and 
Wn = 0.65 (right). The plots cover all the bulk between the two branes. 


[12] Note that the approximation which makes use of an ef¬ 
fective Lagrangian is only compatible with second order 
calculations. 

[13] Note that this conservation equation can be obtained by 
taking the limit y ip in the equation Go4 = 0. 

[14] The unperturbed brane distance is taken to be finite and 
equal to bo/2. 


[15] 

[16] 


In terms of the fundamental scale M one has m = M, 
= M-®, bo = and = M'^. 

Let us recall that this is achieved by adopting the par¬ 
ticular gauge choice for and which sets f^\ (t) = 
fu} it) = 0 and = 0. 
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FIG. 4: The two graphs in the first row show the ratio between to O(e^) and ^ d/a to 0(e) (dotted line) and to O(e^) 
(solid line) at a given time, for p„ = 10~^ and w„ = 0.95 (left) and for p„ = 10“^ AT* and w„ = 0.95 (right). In the two 
graphs in the second row, R^ to O(e^) (solid line) is compared to the modula of to 0(e) (dashed line) and to O(e^) (dotted 
line) for = 2 • 10“^ and Wn = 0.95 (left), and for pn = 6 ■ lO”*^ M'^ and = 0.65 (right). The plots cover all the bulk 
between the two branes. 











